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Borohydride  oxidation  electrokinetics  over  the  Au(l  1 1)  surface  are  simulated  using  first-principles 
determined  elementary  rate  constants  and  a  microkinetic  model.  A  method  to  approximate  the  poten¬ 
tial  dependent  elementary  step  activation  barriers  based  on  density  functional  theory  calculations  is 
developed  and  applied  to  the  minimum  energy  path  for  borohydride  oxidation.  Activation  barriers 
of  the  equivalent  non-electrochemical  reactions  are  calculated  and  made  potential  dependent  using 
the  Butler-Volmer  equation.  The  kinetic  controlled  region  of  the  borohydride  oxidation  reaction  lin¬ 
ear  sweep  voltammogram  over  the  Au(l  1 1)  surface  is  simulated.  The  simulation  results  suggest  that 
B-H  bond  containing  species  are  stable  surface  intermediates  at  potentials  where  an  oxidation  current  is 
observed.  The  predicted  rate  is  most  sensitive  to  the  symmetry  factor  and  the  BH2OH  dissociation  barrier. 
Surface-enhanced  Raman  spectroscopy  confirms  the  presence  of  BH3  as  a  stable  intermediate. 
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1.  Introduction 

Increasing  demand  for  mobile  power  conversion  devices  has 
motivated  research  directed  towards  developing  a  range  of  fuel 
cell  technologies.  Hydrogen  fuel  cells  have  attracted  the  most 
interest  for  large-scale  applications  like  automobiles,  however 
they  are  less  interesting  for  small-scale  power  applications  due 
to  hydrogen  storage  and  gas  handling  safety  concerns.  Alterna¬ 
tively,  direct  borohydride  fuel  cells  (DBFCs)  have  the  potential  to 
convert  chemical  energy  to  electrical  energy  from  a  high  energy 
density  sodium  borohydride  aqueous  fuel,  making  them  attractive 
for  small-scale,  portable  power  applications.  Poor  anode  efficiency 
limits  DBFC  application  [1-4].  The  mechanism  of  direct  oxidation 
and  competing  hydrolysis  reactions  in  alkaline  media  is  unclear 
and  difficult  to  study  experimentally,  therefore,  kinetic  model¬ 
ing  of  the  elementary  reactions  of  borohydride  oxidation  can  help 
to  determine  limiting  steps  and  reactive  intermediates  to  guide 
rational  design  of  improved  catalysts.  We  use  density  functional 
theory  (DFT)  studies  combined  with  microkinetic  modeling  to  pre¬ 
dict  the  potential  dependent  rate  of  borohydride  oxidation  over 
the  Au(l  1 1 )  electrode  surface.  Comparison  between  the  simulated 
and  experimental  linear  sweep  voltammograms  as  well  as  Raman 
spectroscopy  of  adsorbed  intermediates  is  used  to  corroborate  the 
DFT-based  reaction  mechanism. 
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We  previously  studied  the  mechanism  of  borohydride  oxida¬ 
tion  and  hydrolysis  over  gold  and  platinum  surfaces  using  DFT 
methods  [5,6].  We  determined  the  preferred  reaction  path  for 
borohydride  oxidation  over  Au(l  1 1)  and  Pt(  111)  surfaces.  We 
also  identified  the  key  energetic  steps  in  direct  oxidation  and 
competing  hydrolysis  reactions.  Our  previous  analysis  utilized 
DFT-determined  elementary  reaction  thermodynamics,  however, 
reaction  kinetics  were  not  directly  evaluated  due  to  limitations 
in  calculating  potential  dependent  activation  barriers  with  DFT. 
This  motivates  us  to  develop  a  simple  method  to  calculate  the 
potential  dependent  activation  barriers  for  the  borohydride  oxi¬ 
dation  reaction.  Evaluation  of  activation  barriers  will  help  to  clarify 
limiting  elementary  steps  and  determine  energetic  parameters 
which  can  be  evaluated  in  designing  a  new  catalyst  with  better 
activity  and  selectivity  than  gold  and  platinum  electrodes.  With 
approximate  elementary  step  reaction  barriers  from  DFT,  potential 
dependent  elementary  rate  constants  are  calculated  and  input  to 
microkinetic  models.  Microkinetic  models  are  then  used  to  simu¬ 
late  electrokinetic  experiments,  such  as  linear  sweep  voltammetry, 
for  comparison  with  experimental  studies.  DFT  calculations  and 
microkinetic  modeling  may  then  be  coupled  to  predict  the  kinetics 
of  yet  untested  catalyst  formulations,  providing  a  computational 
design  tool. 

For  the  complex  8e_  borohydride  oxidation  reaction,  agreement 
between  DFT-predicted  and  experimental  voltammograms  is  not 
sufficient  to  verify  the  theoretically  determined  reaction  mecha¬ 
nism.  In  this  study,  we  used  surface-enhanced  Raman  spectroscopy 
(SERS)  [7]  detection  of  stable  species  on  roughened  gold  electrodes 
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to  further  validate  the  DFT-predicted  reaction  mechanism.  Gold  is 
not  an  active  catalyst  for  breaking  B-H  bonds  and  we  use  SERS  to 
examine  whether  stable  boron  species  with  SERS-active  B-H  bonds 
are  present  at  various  electrode  potentials. 

In  addition  to  the  specific  borohydride  oxidation  kinetics  stud¬ 
ied  herein,  a  new  method  for  determining  potential  dependent 
activation  barriers  from  DFT  calculations  is  presented.  Atomistic 
modeling  of  elementary  electrocatalytic  rate  constants  is  chal¬ 
lenging,  and  few  methods  exist  to  compute  potential  dependent 
redox  activation  barriers  from  quantum  mechanical  methods.  The 
reaction  center  model  of  Anderson  has  been  applied  to  compute 
potential  dependent  barriers,  however,  it  is  restricted  to  an  unre¬ 
alistic  representation  of  the  electrode  surface  using  a  few  metal 
atoms  [7].  The  double  reference  method  developed  by  Filhol  and 
Neurock  [8,9]  has  been  applied  to  non-redox  elementary  steps  at 
the  electrode  surface  [10],  however,  direct  evaluation  of  barriers 
to  electron  transfer  steps  is  restricted  (currently  applied  in  a  single 
study  [11])  due  to  the  static  nature  of  the  electrolyte  and  limited 
unit  cell  size. 

DFT  evaluation  of  activation  barriers  for  surface  reactions  not 
involving  electron  and  ion  transfer  is  well  established.  Herein,  we 
present  an  approach  to  extrapolate  potential  dependent  barriers  for 
non-electrochemical  hydrogenation/dehydrogenation  reactions  to 
their  equivalent  electrochemical  redox  reactions.  Activation  bar¬ 
riers  for  the  non-electrochemical  reaction  are  assigned  to  an 
appropriate  equilibrium  potential  for  the  elementary  reaction  step. 
The  barriers  are  then  extrapolated  utilizing  the  Butler-Volmer 
assumption  that  the  activation  barrier  of  a  single  redox  reaction 
is  linearly  dependent  on  the  potential  with  a  symmetry  factor  ft. 
For  a  single  oxidation  step: 

Ea(V)  =  E°a-P(V-V0 )  (1) 

where  Ea°  is  activation  barrier  at  the  equilibrium  potential  V0>  and  (3 
is  the  symmetry  factor.  For  most  electrochemical  systems  the  value 
of  p  is  between  0.3  and  0.7  [12]. 

In  this  work,  we  predict  the  rate  of  borohydride  electro¬ 
oxidation  on  the  Au(l  1 1)  surface.  Most  of  previous  experimental 
studies  have  concluded  that  borohydride  oxidation  on  the  gold 
surface  is  a  7-8  electron  process  [13-15],  however,  some  recent 
studies  suggest  that  the  oxidation  might  be  a  four  electron  process, 
also  producing  two  hydrogen  molecules  at  low  potentials  [16,17]. 
These  studies  did  concur  that  the  process  at  higher  potentials  is 
nearly  a  complete  8-electron  process  [16,17].  As  a  first  approxi¬ 
mation,  we  assume  that  borohydride  oxidation  on  the  gold  surface 
is  an  8-electron  process,  and  any  minor  hydrogen  production  at 
lower  potentials  is  neglected.  Though  Au  electrodes  are  reasonably 
selective  to  the  direct  borohydride  oxidation  reaction,  they  suffer 
from  large  overpotentials  due  to  slow  oxidation  kinetics.  The  limit¬ 
ing  elementary  reaction  steps  and  stable  surface  intermediates  are 
determined  based  on  DFT  and  microkinetic  modeling.  A  simulated 
linear  sweep  voltammogram  (LSV)  is  reported  and  compared  to  the 
experimental  LSV.  SERS  results  confirm  stable  intermediates  for  the 
borohydride  oxidation  process  over  the  Au  electrode. 

2.  Research  methods 

2.1.  Experimental  methods 

All  chemicals  were  purchased  from  Sigma-Aldrich  and  used 
as  received  unless  otherwise  stated.  Linear  sweep  voltammetry 
(LSV)  was  performed  using  a  computer-controlled  PARSTAT  2273 
potentiostat  (Princeton  Applied  Research).  A  three-electrode  glass 
cell  (125  mL,  Pine  Instruments)  was  outfitted  with  a  3  mm  gold 
rotating  disk  working  electrode  (Metrohm)  and  a  graphite  rod 
(Sigma-Aldrich)  counter  electrode.  The  working  electrode  was  pol¬ 


ished  to  a  mirror  finish  using  1  pan  and  0.05  pan  alumina  paste 
(Pine  Instruments).  A  single-junction  Ag/AgCl  (4  M  KC1,  Pine  Instru¬ 
ments)  reference  electrode  was  employed  and  2  M  NaOH  served  as 
the  electrolyte.  A  borohydride  concentration  of  0.03  M  was  used 
for  LSV  and  SERS  experiments.  All  the  potentials  reported  here  are 
converted  to  normal  hydrogen  electrode  (NHE)  scale. 

For  SERS  experiments,  the  gold  surface  was  roughened  by 
cycling  for  25  scans  in  a  deoxygenated  aqueous  0.1  M  KC1  solu¬ 
tion  between  -0.08  V  and  1.42  V  with  scan  rate  of  500  mV  s-1 
[18,19].  The  electrode  was  then  cleaned  ultrasonically  for  15  min 
in  18.5  MC2  water.  A  confocal  Raman  microscope  (Renishaw  inVia) 
with  a  20x0.40  NA  glass  objective,  1200  lines  cm-1  diffraction 
grating,  a  solid  state  785  nm  diode  laser,  and  Pelletier  cooled  CCD 
detector  was  used  to  collect  the  SERS  spectra. 

The  working  electrode  was  placed  in  0.03  M  NaBH4  and  2M 
NaOH  solution  at  a  constant  potential  for  1 0  min.  The  electrode  was 
rotated  at  500  rpm  to  avoid  accumulation  of  hydrogen  gas  on  the 
surface.  After  the  1 0  min  potential  hold,  the  electrode  was  removed 
from  solution  under  potential  and  the  SERS  was  collected  imme¬ 
diately  from  the  surface  of  the  electrode.  Co-additions  of  5  scans 
(1950-2250  cm-1)  with  an  exposure  of  10  s  were  used  to  increase 
the  signal  to  noise  ratios  of  the  spectra. 

2.2.  Computational  methods 

Density  functional  theory  calculations  were  performed  using 
the  ab  initio  total  energy  and  molecular-dynamics  Vienna  ab  initio 
simulation  program  (VASP)  developed  at  the  Institute  for  Mate¬ 
rial  Physics  at  the  University  of  Vienna  [20-22].  To  represent  the 
electronic  structure,  the  projector  augmented  wave  method  [23,24] 
was  used  with  a  plane  wave  basis  set  (cut-off  energy  400  eV)  and 
the  Perdew-Wang  [25]  form  of  the  generalized  gradient  approxi¬ 
mation.  To  optimize  the  structures,  a  3  *  3  *  1  Monkhorst-Pack  grid 
was  used  followed  by  a  single-point  calculation  using  a  4  *  4  *  1  grid. 
To  optimize  the  structures,  forces  on  all  atoms  were  minimized  to 
values  lower  than  0.05  eVA-1.  The  climbing  image  nudged  elastic 
band  method  (CI-NEB)  was  used  to  isolate  the  transition  states  of 
elementary  reactions  on  the  surface  [26,27].  Four  images  with  equal 
spacing  were  used  between  the  initial  and  final  state.  The  transition 
state  was  identified  as  the  image  with  the  highest  energy  whereas 
the  tangent  force  is  less  than  0.04  eV  A-1.  Transition  states  were 
confirmed  to  have  a  single  imaginary  vibrational  frequency. 

Similar  to  our  previous  studies  on  borohydride  oxidation  on 
Au(l  11)  [5]  and  Pt(l  11)  [6]  surfaces,  we  used  a  vacuum-slab 
model  of  the  surface.  Ion  adsorption  and  surface  reaction  ener¬ 
gies  as  well  as  optimized  adsorbate  structures  were  previously 
reported  [5].  Our  initial  study  of  boro  hydride  oxidation  on  Au(l  1 1) 
neglected  the  possible  formation  of  a  B(OH)4~  product  or  reac¬ 
tion  through  a  B(OH)3*  surface  intermediate  (where  “*”  is  used 
throughout  to  note  surface  bound  species).  The  binding  energy  (BE) 
ofB(OH)3*  was  calculated  using  the  following  equation: 

B(OH)3+*  ->  B(OH)3* 

BE  =  Edft,B(0H)*  -  £*  -  £dft,b(oh)3  (2) 

where  Edft,b(oh)*  is  the  energy  of  adsorbed  species,  E*  is  the  energy 
of  bare  surface  and  Edft,b(oh)3  is  the  energy  of  isolated  B(OH)3 
species. 

2.2.1.  Calculation  of  potential  dependent  adsorption  equilibrium 
constants 

The  adsorption  of  BH4_  ions  on  the  surface  was  considered  to  be 
at  equilibrium  with  the  bulk  solution  at  all  potentials  of  interest: 

BH4-(aq)  +  H20*  **  BH4*  +  H20(aq)  +  e- 
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BH4-(aq)^  BH4*  (8J  h  BH3*(62)^  BH2*(03) 

K-1  -2  k_3 

^4  ^5  kg  k7 

BH2OH*  (04)  <-►  BHOH*(05)  BH(OH)2*  (06)  <r+  B(OH)2*  (07)  «-* 

k^  k5  k.6  k.7 

B(OH)3*(08)  ^B(OH)4-(aq) 

Fig.  1.  Minimum  energy  path  for  borohydride  oxidation  over  Au(l  1 1)  surface.  Each  reaction  step  produces  an  electron  (other  than  the  final  product  desorption  step)  and 
consumes  a  OH  species  either  to  form  an  H20  molecule  or  a  B-OH  bond. 


For  the  BH4-  concentration  of  0.03  M,  the  adsorption  of  BH4-,  with 
one  molecule  over  nine  Au(l  1 1 )  surface  atoms  is  favorable  (Gibb’s 
free  energy  change  exergonic)  at  potentials  greater  than  0.05VNHE, 
as  determined  in  our  previous  work  [5].  The  Gibbs  free  energy 
change  of  borohydride  adsorption  (AGi),  relative  to  the  0.03 M 
aqueous  concentration,  is  therefore  expressed  as: 

AG,  =(0.05-VNHE)e-  (3) 


where  e~  represents  the  electron  charge  and  is  equal  to  1  when 
using  an  energy  unit  of  eV.  The  equilibrium  constant  for  this  adsorp¬ 
tion  step  (/G[ )  is  then 


K  i 


^BH* 

0h2o*[BH4] 


(4) 


where  0BH*  is  the  surface  fraction  covered  with  BH4*  species  and 
6>H2o*  is  the  surface  coverage  of  water  (replacing  the  typical  free 
surface  coverage  due  to  the  aqueous  environment).  The  concen¬ 
tration  of  the  borohydride  anion  is  included  in  the  denominator, 
though  rigorously  this  should  be  noted  as  the  species  activity.  This 
assumes  an  activity  coefficient  of  1,  and  the  concentration  in  the 
denominator  is  non-dimensionalized  by  the  0.03  M  concentration 
used  in  calculation  of  AGt. 

The  equilibrium  constant  for  the  BH4-  ion  adsorption  is  there¬ 
fore  potential  dependent,  and  we  approximate  that  the  adsorption 
free  energy  is  coverage  independent. 

Similarly,  for  the  adsorption  of  OH-  species  ([OH-]  =  2M,  cal¬ 
culated  for  one  OH  group  per  9  surface  atoms)  the  equilibrium 
constant  is  [5]: 


I<2  =  exp 


0QH* 

£h2o*[OH-] 


(5) 


where 


AG2  =  (0.29  -  Vnhe)^~  (6) 

The  dimensionless  hydroxyl  concentration,  non-dimensionalized 
by  the  reference  concentration  of  2  M,  is  included  in  the  denomi¬ 
nator  of  Eq.  (5)  rather  than  the  more  rigorous  activity  expression. 


2.2.2.  Calculation  of  the  desorption  equilibrium  constant  of  boric 
acid 

We  have  previously  concluded  that  boric  acid  is  the  final  prod¬ 
uct  of  the  surface  borohydride  oxidation  reaction  over  the  Au(  111) 
electrode  [6].  The  product  boric  acid  may  then  reaction  in  the  solu¬ 
tion  phase  to  produce  the  hydrated  borate  anion: 


the  free  energy  change  of  replacing  adsorbed  B(OH)3*  with  H20 
followed  by  the  equilibrium  of  Eq.  (7): 


AG3  =  Gh20*  +  GB(OH)-  “  GOH-q  “  GH2Oaq  “  GB(OH);  (8) 

The  calculated  value  of  AG3  is  -0.74  eV. 

The  evaluation  of  the  free  energy  of  the  solvated  and  surface 
bound  species  is  discussed  in  our  previous  publications  [5,6]. 

The  equilibrium  constant  for  the  B(OH)3*  desorption  step  is 
then: 


I<3  =  exp 


^h2o*[B(OH)4  ] 

^B(OH)3*[°H  1 


(9) 


2.2.3.  Potential  dependent  activation  barriers 

Density  functional  theory  calculations  do  not  trivially  pro¬ 
vide  potential  dependent  activation  barriers  for  redox  reactions. 
We  calculated  non-potential  dependent  activation  barriers  for  the 
equivalent  chemical  reactions  on  the  Au(l  1 1 )  surface,  and  extrap¬ 
olated  these  to  be  potential  dependent  using  a  Butler-Volmer 
formalism.  The  overall  oxidation  reaction  is  presumed  to  proceed 
through  the  minimum  energy  path,  and  competing  less  favorable 
elementary  steps  are  not  considered. 

The  elementary  oxidation  steps  on  the  Au(l  1 1)  surface  could 
be  categorized  to  three  different  types  of  reaction:  (1)  B-H  cleav¬ 
age  reactions,  (2)  B-0  bond  formation  reactions  and  (3)  O-H  bond 
cleavage  reactions.  Breaking  O-H  bonds  is  not  favorable  as  gold  is 
a  poor  catalyst  for  activating  these  bonds,  and  we  conclude  that 
the  final  eight  electron  oxidation  product  of  the  surface  reaction  is 
boric  acid,  B(OH)3.  Alternatively,  B-0  formation  reactions  are  very 
favorable  and  exothermic  on  the  gold  surface  and  our  DFT  calcu¬ 
lations  indicate  that  the  activation  barriers  of  these  reactions  are 
negligible.  Activating  B-H  bonds  on  the  surface  limits  the  overall 
oxidation  process. 

Fig.  1  shows  the  minimum  energy  path  for  borohydride  oxida¬ 
tion  on  the  Au(  111)  surface.  In  Fig.  1 ,  kx  represents  the  rate  constant 
for  each  elementary  electrooxidation  reaction  and  /<_*  represents 
the  reverse  reduction  rate  constant.  The  OH-,  H20,  and  electron 
species  involved  in  each  reaction  step  are  left  off  of  Fig.  1  for  sim¬ 
plicity.  Each  reaction  step,  other  than  the  final  desorption  step, 
produces  an  electron.  All  reaction  steps,  other  than  the  initial  BH4- 
adsorption,  consume  a  OH-  species  to  either  form  a  water  molecule 
or  a  B-OH  bond. 

For  a  typical  B-H  breaking  step  on  the  surface,  the  non¬ 
electrochemical  dehydrogenation  elementary  reaction  can  be 
written  as: 

A-H*  ->►  A*  +  H*  (10) 


B(OH)3(aq)  +  OH-(aq)  B(OH)4-(aq)  (7) 

The  surface  adsorbed  boric  acid  species,  B(OH)3*,  is  assumed  to  be 
in  equilibrium  with  the  B(OH)4-  ion  in  the  solution  phase.  The  free 
energy  of  desorption  of  B(OH)3*  may  be  calculated  by  evaluating 


The  equivalent  electro-oxidation  step  for  this  reaction  is: 

A-H*  +  OH-aq  ->  A*  +  H2Oaq  +  e-  (11) 

For  the  electrocatalytic  reaction,  the  B-H  bond  is  activated  through 
interaction  with  the  catalyst  surface.  We  approximate  that  the 


G.  Rostamikia  et  al.  /  Journal  of  Power  Sources  196(2011)  9228-9237 


9231 


activation  barriers  of  these  two  reactions  are  equivalent  at  the  elec¬ 
trochemical  potential  at  which  the  chemical  potential  of  the  H* 
product  is  in  equilibrium  with  its  desorbed  product: 

H*  +  OH-aq  H2Oaq  +  e-  +  *  (12) 


*1010 

[bh4-] 


09 


*2  0io 

[OH-] 


(24) 

(25) 


With  this  assumption,  the  activation  barrier  for  the  non¬ 
electrochemical  dehydrogenation  reaction  is  assigned  to  the 
potential  at  which  the  surface  H*  species  has  the  same  chemical 
potential  as  the  redox  reaction  product.  At  this  potential,  the  dis¬ 
sociated  hydrogen  atom  has  equal  thermodynamic  preference  to 
adsorb  to  the  surface  or  oxidize  to  produce  a  proton  and  an  electron. 

The  activation  barriers  for  the  non-electrochemical  reaction  10 
are  calculated  with  the  vacuum-slab  DFT  model.  Then,  similar  to 
the  Butler-Volmer  expression,  potential  dependent  activation  bar¬ 
riers  for  oxidation  reactions  are  linearly  related  to  the  equivalent 
electrochemical  reaction  activation  barriers  using  a  symmetry  fac¬ 
tor.  For  example,  the  activation  barrier  of  oxidative  dissociation  of 
BH4*  (Efl( V)),  BH4*  +  OH"(aq)  BH3*  +  H20(aq)  +  e 
is  written  as: 

Ea,(V)  =  Eai°-P(V-V 0)  (13) 

where  Efll  0  is  the  activation  barrier  for  the  chemical  step,  f3  is  the 
symmetry  factor,  and  Vo  is  the  equilibrium  potential  for  the  oxida¬ 
tion  of  surface  hydrogen. 


2.2.4.  Calculation  of  potential  dependent  rate  constants 

Using  transition  state  theory,  an  Arrhenius  form  is  considered 
for  the  elementary  rate  constants.  For  example,  for  the  initial  B-H 
dissociation  step: 

ki=Aj*e~Ea  '(V)/RT  (14) 


The  pre-exponential  factor  is  approximated  as: 


(15) 


where  kb  is  the  Boltzman  constant,  T  is  the  absolute  temperature, 
h  is  Planck’s  constant,  and  r  represents  the  probability  that  the 
system  will  go  over  the  barrier.  A  reasonable  value  of  r  =  0.5  is  used 
in  all  calculations.  The  rate  of  any  surface  reaction  is  calculated  as: 


-rn=l<n0n  (16) 

where  rn  represents  the  rate  of  nth  reaction,  kn  is  the  rate  constant 
for  the  nth  step,  and  0n  is  the  surface  fraction  of  the  species  being 
converted  in  step  n.  The  units  of  the  rate  calculated  are  s-1  site-1. 
Since  each  elementary  reaction  produces  1  electron,  the  current 
will  be  in  units  of  electrons  s-1  site-1. 


2.2.5.  Micro-kinetic  model  construction 

A  set  of  algebraic  and  differential  equations  are  written  as: 


^  =  fc101-k202  +  k_203-k_i02 

(17) 

df) 

=  k202  -  I0O3  -  k_: 203  +  k.304 

(18) 

=  k3#3  —  k404  —  k_3$4  +  k- 4$5 

(19) 

^  =  k4e4  -  k5e 5  -  k_4e5  +  k_  5e6 

(20) 

df) 

-gjT  =  &505  -  ^606  -  k-506  +  1<_q0j 

(21) 

—^r  =  1<q0q  -  l<j0j  -  k_607  +  k_7$8 

(22) 

M 

II 

(23) 

where  0\  to  08  were  introduced  in  Fig.  1  and  correspond  to  the  sur¬ 
face  fraction  of  different  boron  containing  intermediate  species.  69 
is  the  surface  fraction  of  OH  species  and  6\q  is  the  “empty”  sur¬ 
face  fraction,  or  better  written  in  solution  as  the  fraction  of  surface 
which  is  covered  with  water.  Since  the  concentration  of  B(OH)4- 
ions  is  very  small,  lxl  0-6  M  is  assigned  for  the  B(OH)4-  concentra¬ 
tion,  and  this  is  used  to  evaluate  the  B(OH)3*  surface  concentration 
as  described  above.  The  value  of  1  x  10-6  M  represents  an  approx¬ 
imate  concentration  value  for  the  B(OH)4-  ion  given  mass  transfer 
into  the  stagnant  solution  from  the  surface  region  [28].  Though  this 
value  is  relatively  arbitrary,  the  exact  value  has  little  impact  on 
the  simulated  kinetics  as  product  adsorption  is  negligible  and  the 
reverse  reduction  reaction  rates  are  miniscule  at  the  high  oxidation 
overpotentials  necessary  to  observe  an  anodic  current. 

In  the  set  of  differential  equations,  all  rate  constants  and  equi¬ 
librium  constants  are  dependent  on  the  electrode  potential.  The 
electrode  potential  can  then  be  made  time  dependent  given  an 
initial  potential  and  sweep  rate  in  the  LSV  experiment.  The  set  of 
algebraic  and  differential  equations  can  be  solved  simultaneously 
to  find  the  time/potential  dependent  current  of  borohydride  oxida¬ 
tion  over  the  Au(l  1 1 )  surface.  Numerical  solution  requires  a  small 
time  step  to  reach  reasonable  real  times  because  both  extremely 
fast  and  slow  reactions  are  included  in  the  model.  As  a  first  approx¬ 
imation,  we  assume  that  the  surface  reactions  reach  a  quasi-steady 
state  as  the  potential  is  varied.  Because  the  time  derivative  of  sur¬ 
face  intermediates  is  approximated  as  zero,  dOnldt  is  equal  to  zero. 
This  leaves  a  set  of  algebraic  equations  which  can  be  solved  at 
each  potential.  The  potential  dependent  fractional  coverages  of  all 
species  (0n)  are  determined  and  used  to  calculate  the  reaction  rate, 
as  described  in  the  next  section. 


2.2.6.  Calculating  the  current  from  the  microkinetic  model 

As  the  overall  oxidation  process  produces  eight  electrons,  the 
overall  current  could  be  calculated  by  multiplying  the  rate  of  any 
surface  oxidation  reaction  by  eight.  For  example: 

ic  =  8 k7q7  in  units  of  electrons  s-1  surface  site-1 


To  convert  the  per  site  current  to  a  per  area  current,  we  multi¬ 
ply  by  the  electronic  charge  and  the  number  of  surface  sites  per 
area.  The  electrode  used  experimentally  has  a  surface  area  (A)  of 
7.065  x  10-2  cm2.  For  the  Au(l  1 1)  surface,  the  area  occupied  by 
one  atom  (a)  is  7.0298  x  10-16  cm2.  The  number  of  atoms  on  the 
electrode  surface  is: 


(26) 


If  we  assume  that  n  atoms  make  one  reaction  site  on  the  surface, 
then  the  number  of  active  sites  (S)  is: 


(27) 


A  value  of  n  =  3  is  assigned  for  number  of  surface  atoms  per  active 
site.  This  value  is  based  on  the  size  of  intermediates,  as  all  surface 
intermediates  are  similar  in  size  to  three  Au  surface  atoms.  The 
effect  of  changing  n  is  discussed  in  the  results  section. 

The  computationally  calculated  current  (/)  is  given  by: 


/  =  8 eSk797  (  j) 


where  e  is  the  electronic  charge. 


(28) 
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Fig.  2.  Adsorption  configuration  of  B(OH)3  on  the  Au(l  1 1)  surface. 


3.  Results  and  discussion 

The  electro-oxidation  of  borohydride  is  modeled  over  the 
Au(l  1 1)  surface.  As  our  previous  paper  applying  DFT  to  the  study 
of  borohydride  oxidation  over  the  Au(l  1 1)  surface  neglected  the 
possibility  of  forming  boric  acid  the  final  8e~  product,  the  preferred 
binding  site  of  boric  acid  is  first  reported  in  Section  3.1.  The  acti¬ 
vation  barriers  for  the  favorable  oxidation  path  are  presented  in 
Section  3.2.  The  simulated  borohydride  oxidation  LSV  results  are 
presented  in  Section  3.6  and  compared  with  experiment.  Raman 
spectroscopy  results  are  presented  in  Section  3.7  for  further  cor¬ 
roboration  of  the  computationally  predicted  reaction  mechanism. 


3.1.  Preferred  binding  site  of  boric  acid  on  the  Au(  111)  surface 

The  binding  energy  and  optimal  binding  configuration  of  boro¬ 
hydride  oxidation  intermediates  was  previously  reported  [5].  This 
study  did  not  consider  production  of  a  final  boric  acid  (B(OH)3)* 
product,  which  has  since  been  determined  as  favored  over  B02* 
as  the  final  product  of  the  surface  oxidation  reaction.  Fig.  2  illus¬ 
trates  the  most  stable  adsorption  configuration  for  B(OH)3*  on  the 
Au(l  1 1)  surface.  B(OFI)3*  species  are  stable  molecules  and  have  a 
weak  interaction  with  the  Au(l  1 1)  surface.  The  binding  energy 
of  these  species  to  the  Au(l  1 1)  surface  is  -0.3  eV.  The  desorp¬ 
tion  of  these  species  in  aqueous  media  therefore  has  only  a  slight 
activation  barrier  and  will  likely  be  mass  transfer  limited. 


3.2.  Non-electrochemical  activation  barriers  for  the  preferred 
reaction  path 

The  activation  barrier  for  the  initial  dehydrogenation  reaction 
was  previously  reported  [5]: 

BH;  BH5  +  H*  Eact  =  0.37  eV  (36  kj  mol-1 )  (29) 

The  activation  barriers  for  the  B-0  formation  reactions 
BH2*  +  OH*^  bh2oh* 

BHOH*  +  OH*^  BH(OH)2* 

B(OH)2*  +  OH*^  B(OH)3* 

were  evaluated  and  all  found  to  be  extremely  small  (<0.01  eV).  For 
microkinetic  modeling,  a  barrier  of  0.01  eV  is  assigned  for  all  B-0 
formation  steps. 

The  transition  state  for  BFI3*  dehydrogenation  to  BH2*  and 
FI*  is  illustrated  in  Fig.  3.  The  barrier  for  this  reaction  is  0.81  eV 
(78kJmol_1).  The  dissociating  B-FI  bond  stretches  from  1.23  A  in 
the  BH3*  state  to  2.3  A  at  the  transition  state.  A  single  imaginary 
vibrational  frequency  of  i284  cm-1  confirms  that  the  state  in  Fig.  3b 
is  a  transition  state. 

The  equilibrium  and  transition  states  for  BFI2OFI*  dehydrogena¬ 
tion  to  BFIOFI*  and  FI*  are  illustrated  in  Fig.  4.  This  step  has  the 
highest  activation  barrier  among  the  B-FI  dissociation  steps.  The 
activation  barrier  is  1 .06  eV  ( 1 02  kj  mol-1 ).  The  B-H  bond  extends 
from  1.2  A  at  the  initial  state  to  1.9  A  at  the  transition  state.  An 
imaginary  vibrational  frequency  of  i'406  cm-1  confirms  that  Fig.  4b 
represents  a  transition  state. 

The  transition  state  for  BH(OH)2*  dissociation  to  B(OH)2*  and  H* 
is  illustrated  in  Fig.  5.  Similar  to  BFI2OFI*  dissociation,  this  step  has  a 
high  barrier  (0.91  eV  or  88  kj  mol-1 ).  The  dissociating  bond  extends 
from  1 .2  A  at  the  initial  state  to  2.3  A  at  the  transition  state.  An  imag¬ 
inary  vibrational  frequency  of  i282  cm-1  confirms  the  transition 
state. 

Table  1  lists  the  activation  barriers  for  B-FI  dehydrogenation  and 
B-0  formation  reactions  on  the  Au(l  1 1 )  surface. 

3.3.  Surface  hydrogen  oxidation  potential  on  the  Au(  111)  surface 

The  activation  barriers  calculated  for  dehydrogenation  reactions 
are  assigned  to  a  specific  potential  by  approximating  that  they  rep¬ 
resent  the  appropriate  barrier  at  the  chemical  potential  for  which 
the  product  FI*  species  is  in  equilibrium  with  solution,  as  expressed 
in  Eq.  (12).  Assuming  that  Eq.  (12)  is  at  equilibrium,  for  a  1/9 
monolayer  of  hydrogen  atoms  on  the  surface  and  pFI  of  14.26,  the 


Fig.  3.  BH3*  dehydrogenation  to  BH2*  and  IT  over  the  Au(l  1 1)  surface:  (a)  initial  state,  (b)  transition  state,  and  (c)  final  state. 
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Fig.  4.  Dehydrogenation  of  BH2OH*  to  BHOH*  and  H*  over  the  Au(l  1 1)  surface:  (a)  initial  state,  (b)  transition  state,  and  (c)  final  state. 


Fig.  5.  Dehydrogenation  of  BH(OH)2*  to  B(OH)2*  and  H*  over  the  Au(l  1 1)  surface:  (a)  initial  state,  (b)  transition  state,  and  (c)  final  state. 


equilibrium  potential  for  surface  hydrogen  oxidation  is  calculated 
via: 

V0(NHE )  =  GH2oaq  +  C*  -  CH*  -  C0H-  -  4.6  (30) 

aq 

where  4.6  is  conversion  factor  between  absolute  (vacuum  scale) 
potential  and  the  normal  hydrogen  electrode  scale  [5]  determined 
using  DFT  methods  consistent  with  those  used  herein.  The  value  of 
V0  calculated  using  Eq.  (30)  is  -1.19  V. 

3.4.  Potential  dependent  rate  constants  for  elementary  oxidation 
reactions 

Given  the  non-electrochemical  barriers  in  Table  1  and  the 
calculated  equilibrium  potential  for  the  H*  desorption,  potential 
dependent  rate  constants  for  elementary  reactions  are  calculated 
via  Eq.  (14).  The  potential  dependent  activation  barriers  for  all  steps 
are  calculated  using  Eq.  (13),  with  V0  equal  to  -1.19V(NHE)  and 
Ea°  from  Table  1.  Table  2  gives  the  potential  dependent  activation 
barriers  for  all  elementary  forward  reaction  steps  of  the  minimum 
energy  path. 


For  the  backward  (reduction)  reactions  the  potential  dependent 
activation  barriers  is  calculated  similarly.  For  example: 

Ealh(V)=  Ea\b+oc{V  +  1.19) 

where  a  =  1  -  is  the  symmetry  factor  of  the  backward  reaction 

[12]. 

The  pre-exponential  factor  for  the  rate  constants  is  calcu¬ 
lated  via  Eq.  (15).  At  room  temperature  (298  K),  the  calcu¬ 
lated  pre-exponential  factor  is  3.1  x  1012  s-1  site-1.  An  average 
value  of  r  equal  to  0.5  is  used  within  the  pre-exponential 
factor. 

Fig.  6  presents  the  elementary  reaction  energetics  and  activation 
barriers  of  elementary  reactions  for  the  preferred  oxidation  path  at 
-0.2  V  and  -IV.  A  value  of  /3  =  0.44  is  considered  for  symmetry 
factor  and  is  discussed  in  Section  3.6.  Fig.  6  shows  that  at  -1 V,  ini¬ 
tial  BH4-  adsorption  is  not  favorable  and  B-H  bond  dissociation 
steps  present  high  barriers.  At  -0.2  V,  all  barriers  become  acces¬ 
sible  at  room  temperature  and  the  overall  oxidation  reaction  will 
proceed. 


Table  1  Table  2 

Activation  barriers  (eV)  for  elementary  B-H  dissociation  and  B-0  formation  reac-  Potential  dependent  activation  barriers  (eV)  for  elementary  oxidation  reactions  of 
tions  and  the  backward  reactions  over  the  Au(l  1 1 )  surface.  borohydride  over  the  Au(l  1 1 )  surface. 


Reaction 

Ea  forward  reaction 

Ea  backward  reaction 

Reactions 

Ea  (V) 

(1)BH4*«>BH3*  +  H* 

£a  if  =  0.37 

Ea  ib  =  0.16 

( 1 )  BH4*  +  OH-  BH3*  +  H20  +  e- 

£a if  =  0.37  -  /3  (V  + 1.19) 

(2)  BH3*  **  BH2*  +  H* 

Ea2f  =  0.81 

Ea2b  =  0.48 

(2)  BH3*  +  OH-  BH2*  +  H20  +  e- 

£a2f  =  0.81  -  /3  (V  + 1.19) 

(3)  BH2*  +  OH*  BH2OH* 

£a3f  =  Negligible  (<0.01 ) 

£a3b  =  2.63 

(3)  BH3*  +  OH-  BH2OH*  +  e- 

£a3f  =  Negligible  (<0.01 ) 

(4)  BH2OH**>BHOH*  +  H* 

£a4f  =  1-06 

£a4b  =  0.67 

(4)  BH2OH*  +  OH”  BHOH*  +  H20  +  e~ 

£a4f=  1.06  -P(V+ 1.19) 

(5)  BHOH*  +  OH*  **  BH(OH)2* 

£a5f  =  Negligible 

£a  5b  =  2.60 

(5)  BHOH*  +  OH-  **  BH(OH)2*  +  e~ 

Ea  5f  =  Negligible 

(6)  BH(OH)2*  B(OH)2*  +  H* 

£a6f=0.91 

Ea  6b  =  0.78 

(6)  BH(OH)2*  +  OH-  B(OH)2*  +  H20  +  e~ 

£a6f  =  0.91  -  /3  (V+1A9) 

(7)  B(OH)2*  +  OH*  **  B(OH)3* 

£a7f  =  Negligible 

Ea7h  =  2.72 

(7)  B(OH)*  +  OH-  B(OH)3*  +  e~ 

Ea  7f  =  Negligible 

9234 


G.  Rostamikia  et  al.  /  Journal  of  Power  Sources  196  (201 1 )  9228-9237 


Fig.  6.  Reaction  energetics  and  activation  barriers  of  the  borohydride  oxidation  preferred  path  at  two  different  potentials,  f3  =  0.44,  [BH4-  ]  =  0.03  M,  [OH-  ]  =  2  M,  T= 298. 


3.5.  Experimental  linear  sweep  voltammetry 

Linear  sweep  voltammetry  data  for  borohydride  oxidation  is 
presented  in  Fig.  7  for  different  scan  rates  (5,  10  and  25mVs-1). 
At  potentials  higher  than  -0.35  V,  the  reaction  rate  is  affected  by 
mass  transfer  and  the  current  varies  with  scan  rate.  At  potentials 
lower  than  -0.35  V,  mass  transfer  has  minimal  effect  on  the  current, 
the  current  is  nearly  independent  of  scan  rate,  and  the  oxidation 
process  is  controlled  by  kinetic  limitations.  Similar  currents  for  dif¬ 
ferent  scan  rates  also  suggest  that  the  surface  reactions  reach  a 
quasi-steady  state  at  each  potential  as  the  potential  is  scanned. 
Therefore,  at  potentials  lower  than  -0.35  V,  we  may  approxi¬ 
mate  that  the  current  is  not  time  dependent  and  only  depends  on 
potential. 


Fig.  7.  Linear  sweep  voltammogram  of  borohydride  oxidation  over  a  gold  electrode, 
[BH4-]  =  0.03 M,  [OH-]  =  2 M,  T=298K,  electrode  diameter  =  3 mm,  scan  rates:  (a) 
25mVs-1,  (b)  10mVs-1,  and  (c)  5mVs-1. 


3.6.  Simulated  linear  sweep  voltammogram 

Eqs.  (19)-(25)  were  solved  at  each  potential  to  determine  the 
simulated  current  as  a  function  of  potential.  Different  values  of 
symmetry  factor  were  considered.  Based  on  experimental  LSV 
results,  the  simulated  current  was  limited  to  values  lower  than 
400  [jlA  to  avoid  consideration  of  transport  limitations  in  the  sys¬ 
tem.  Fig.  8  compares  the  experimental  LSV  data  and  simulation 
results  for  three  different  symmetry  factor  values.  A  single  symme¬ 
try  factor  was  used  for  all  reaction  steps;  the  evaluation  of  symme¬ 
try  factors  directly  from  DFT  calculations  will  be  a  subject  of  future 
study.  Data  in  Fig.  8  assumes  3  surface  atoms  as  one  active  site. 

For  p  values  of  0.5,  the  simulation  overestimates  the  oxida¬ 
tion  current.  The  value  of  0.5  is  a  reasonable  approximate  value 


Fig.  8.  Simulated  (solid)  and  experimental  (dashed)  linear  sweep  voltammo- 
grams  of  borohydride  over  the  gold  electrode.  Different  values  of  symmetry  factor 
are  considered,  [BH4-]  =  0.03 M,  [OH-]  =  2M,  7=298K,  disk  diameter  =  3 mm,  scan 
rate  =  5  mV  s-1 :  (a)  ft  =  0.5,  (b)  ft  =  0.44,  (c)  ft  =  0.4. 
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Fig.  9.  Computational  (solid)  and  experimental  (dashed)  data  for  linear  sweep 
voltammetry  of  borohydride  over  the  gold  electrode.  Effect  of  assuming  different 
values  of  surface  atoms  per  surface  site  in  the  model  is  presented,  [BH4-]  =  0.03  M, 
[OH- ]  =  2  M ,T=  298  K,  disk  diameter  =  3  mm,  scan  rate  =  5  mV  s-1 :  (a)  n  =  2,  (b)  n  =  3, 
(c)  n  =  4. 

in  absence  of  actual  data,  however,  typical  values  fall  between 
0.3  and  0.7  [12].  For  equal  to  0.44,  experimental  and  simulation 
results  are  in  agreement.  Clearly,  the  take-off  potential  and  oxida¬ 
tion  current  are  highly  sensitive  to  the  value,  and  the  need  to 
approximate  this  value  is  a  limitation  of  the  methods  used.  Devi¬ 
ations  from  0.5  can  be  approximated  by  considering  the  variable 
interactions  of  the  equilibrium  and  transition  state  dipole  moments 
with  the  interfacial  electric  field,  and  this  will  be  a  subject  for  future 
studies. 

Fig.  9  presents  the  effect  on  n  (number  of  surface  atoms  per  sur¬ 
face  site)  on  the  computational  current  calculated  from  the  kinetic 
model.  Based  on  the  size  of  the  reaction  intermediate  molecules 
we  initially  considered  a  reasonable  value  of  n  equal  to  3  in  our 
calculation.  Choosing  different  n  values  changes  the  absolute  mag¬ 
nitude  of  the  I-V  curve.  For  a  value  of  n  equal  to  2,  the  I-V  curve 
shifts  to  lower  potentials  due  to  a  higher  number  of  surface  sites 
per  physical  area.  A  value  of  n  equal  to  2  is  likely  not  reasonable 
due  to  the  size  of  the  intermediates  molecules.  Increasing  n  to  4 
decreases  the  current  at  a  given  potential.  The  effect  of  the  n  value 
on  the  computational  current  is  lower  than  the  effect  of  symmetry 
factor  (3,  as  reasonable  n  values  are  not  beyond  the  2-4  range. 

Fig.  10  shows  the  surface  coverage  of  boron  containing  inter¬ 
mediates  {0b  summed  over  all  intermediates)  on  the  Au(l  1 1) 
surface  as  a  function  of  electrode  potential.  The  modeling  results 
suggest  that  BFI2OH*  is  the  most  abundant  surface  intermedi¬ 
ate  and  the  Au(l  1 1)  surface  will  be  nearly  completely  covered 
by  this  species  at  potentials  for  which  the  oxidation  current  is 
substantial.  At  the  potential  of  -0.55  V,  the  surface  fractional 
coverage  of  BH2OH  exceeds  0.99  indicating  essentially  a  full  mono- 
layer  (1  BH2OH  species  per  3  Au  surface  atoms).  At  potentials 
higher  than  -0.55  V  the  oxidation  reaction  is  limited  by  the  oxi¬ 
dation  rate  of  this  species.  This  is  a  direct  consequence  of  the 
BH2OH  oxidation  reaction  barrier  being  larger  than  for  all  other 
elementary  steps. 

The  predicted  rate  of  reaction  is  highly  sensitive  to  the  BH2OH 
dissociation  barrier.  Increasing  the  BH2OH  dissociation  barrier 
decreases  the  predicted  rate  of  reaction;  however,  increasing  or 
decreasing  the  activation  barriers  for  other  elementary  steps  by  5% 
does  not  affect  the  predicted  rate  of  oxidation. 


V(NHE) 


Fig.  10.  Surface  coverage  of  boron  containing  intermediate  on  the  Au(l  1 1 )  surface 
as  a  function  of  potential  from  the  kinetic  model;  n  =  3,  ft  =  0.44. 


Raman  Shift  (cm'1) 

Fig.  11.  Surface  enhanced  Raman  spectra  for  borohydride  oxidation  over  gold  elec¬ 
trode  at  different  potentials. 

3.7.  Surface  enhanced  Raman  spectroscopy 

Fig.  11  shows  surface  enhanced  Raman  spectra  on  a  gold 
electrode  following  borohydride  oxidation  at  various  constant 
potentials  held  for  10  min.  Two  peaks  are  observed  in  the  Raman 
spectra  at  approximately  2125  cm-1  and  2225  cm-1  after  oxidation 
at  -0.4  and  -0.2  V,  which  are  attributed  to  a  stable  intermediate 
adsorbate  on  the  electrode  surface.  No  peaks  are  observed  asso¬ 
ciated  with  borohydride  species  at  the  lower  potentials  of  -0.8  V. 
Because  of  the  negligible  barrier  of  B-0  bond  formation  reactions, 
BH2\  BHOFI*  and  B(OFI)2*,  are  ruled  out  as  possible  stable  sur¬ 
face  intermediates.  Four  surface  species,  BFI4*,  BFI3*,  BH2OFI*  and 
BH(OFI)2*,  require  B-H  dissociation  reactions  for  further  conversion 
and  are  possible  as  stable  intermediates.  The  vibrational  frequen¬ 
cies  of  these  surface  species  were  calculated  using  DFT  methods. 
The  functional  used  in  the  DFT  method  (PW91 )  is  known  to  system¬ 
atically  overestimate  vibrational  frequencies,  and  a  scaling  factor 
of  0.96  was  used  to  correct  the  calculated  frequencies  [29]. 
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Table  3 

Corrected  vibrational  frequencies  (cm-1)  in  the  B-H  stretching  range  for  possible 
stable  intermediates  over  the  Au(l  1 1)  surface.  A  full  list  of  vibrational  frequencies 
of  the  adsorbed  species  is  included  in  the  supplementary  information. 


bh4* 

BH2OH* 

BH(OH)2* 

bh3* 

2335 

2470 

2461 

2231 

1932 

2320 

2204 

1923 

2124 

1868 

The  assignment  of  the  observed  Raman  peaks  to  BH4*,  BH2OH*, 
or  BH(OH)2*  adsorbates  can  be  ruled  out  by  comparison  with 
calculated  frequencies.  The  corrected  B-H  bond  stretching  frequen¬ 
cies  for  BH4*  species  are  2335  cm-1,  1932  cm-1,  1923  cm-1  and 
1868  cm-1  which  do  not  closely  match  the  observed  frequencies 
in  the  Raman  experiments.  Due  to  the  lower  barrier  of  BH4  dis¬ 
sociation  on  the  surface,  these  species  are  not  observed  as  stable 
intermediates  in  the  SERS  experiments.  Similarly,  the  observed 
peaks  in  the  SERS  spectra  could  not  correspond  to  BH2OH  or 
BH(OH)2  species  on  the  surface.  Due  to  the  very  weak  binding 
of  these  stable  species  on  the  Au(l  1 1)  surface  [5],  any  unreacted 
species  of  BH2OH  and  BH(OH)2  could  desorb  from  the  electrode  sur¬ 
face  upon  removing  the  electrode  from  the  reaction  environment. 
Therefore,  these  species  are  not  detected  in  the  SERS  experiments. 

The  corrected  B-H  vibrational  frequency  of  BH3  species  on  the 
Au(l  1 1)  surface  are  2231  cm-1, 2204  cm-1  and  2124  cm-1  and  rea¬ 
sonable  agree  with  the  experimentally  observed  spectra  in  Fig.  11. 
The  substantial  barrier  of  BH3  dissociation  and  relatively  stronger 
binding  (-0.55  eV)  of  these  species  on  the  surface  [5],  allow  detec¬ 
tion  of  the  BH3*  intermediates  as  stable  species  on  the  surface.  The 
observed  Raman  peaks  could  not  be  assigned  to  stable  BH3OH* 
since  corrected  B-H  vibrational  frequencies  for  these  species  are 
all  lower  than  2000  cm-1.  The  loss  the  higher  frequency  peak  at 
-0.1  V  and  broadening  of  the  lower  frequency  peak  has  not  been 
explained,  and  may  represent  a  change  in  coverage  at  higher  oxi¬ 
dation  rates. 

All  vibrational  calculations  considered  intermediates  at  the 
metal-vacuum  interface,  and  the  presence  of  solvating  water  or 
surface  electric  fields  may  alter  the  absolute  values  of  the  vibra¬ 
tional  frequencies.  To  quantify  the  vibrational  shift  induced  due 
to  hydrogen  bonding  with  nearby  water  molecules,  we  calculated 
the  vibration  frequency  of  BH3*  with  3  water  molecules  around 
the  adsorbed  borane  species.  The  corrected  vibration  frequencies 
of  B-H  bonds  at  the  desired  region  changed  from  2231,  2204  and 
2124  cm-1  to  2256, 2249  and  2163  cm-1 ,  respectively  (see  Table  3). 

To  the  best  of  our  knowledge  there  is  only  one  in  situ  infrared 
study  of  the  borohydride  oxidation  reaction  over  the  gold  elec¬ 
trode.  Concha  and  coworkers  reported  a  band  at  1184  cm-1  which 
they  assigned  to  a  B-H  bond  of  either  BH3  or  BH3OH-  species  [30]. 
They  also  reported  two  other  bands  at  1326  and  1415  cm-1  which 
they  assigned  to  the  B-0  bond  of  B02-  or  BOB  species.  We  did 
not  considered  B-O-B  bond  formation  of  BH3OH-  species  in  our 
analysis. 

4.  Conclusions 

A  simple  model  was  developed  to  calculate  the  potential  depen¬ 
dent  activation  barriers  of  borohydride  oxidation  over  the  Au(  111) 
surface.  The  kinetic  region  of  the  linear  sweep  voltammogram 
was  simulated  using  density  functional  theory  calculated  reaction 
energetics  and  agrees  with  the  experimental  voltammogram.  The 
predicted  rate  is  sensitive  to  the  electrochemical  reaction  symme¬ 
try  factor  and  the  BH2OH  dissociation  barrier.  A  value  of  0.44  is 
suggested  for  symmetry  factor  of  borohydride  oxidation  reaction 
over  the  Au(l  1 1 )  surface.  Both  calculated  and  experimental  kinet¬ 
ics  are  consistent  with  a  single  electron  transfer  limiting  reaction 


step,  which  is  determined  to  be  associated  with  B-H  bond  dissoci¬ 
ation  of  a  surface  intermediate. 

Due  to  the  substantial  barrier  of  B-H  dissociation,  BH3*,  BH2OH* 
and  BH(OH)2*  are  possible  stable  surface  intermediates  of  the  reac¬ 
tion,  however,  the  LSV  simulation  suggests  the  surface  is  covered 
by  a  monolayer  of  BH2OH*  at  potentials  higher  than  -0.55 V{NHE). 
SERS  was  used  to  detect  stable  BH3  intermediates  on  the  electrode 
surface  after  oxidation  at  -0.2  and  -OAV(NHE).  Comparison  of  SERS 
spectra,  taken  following  removal  of  the  electrode  from  the  electro¬ 
chemical  cell,  with  DFT-calculated  vibrational  frequencies  suggests 
a  stable  BH3*  intermediate.  The  theory  and  experimental  results 
are  in  qualitative  agreement  that  stable  intermediates  are  formed 
during  borohydride  oxidation  containing  B-H  bonds,  further  sug¬ 
gesting  the  rate  of  borohydride  electrooxidation  is  limited  on  Au 
electrodes  by  slow  B-H  dissociation. 
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Appendix  A.  List  of  variables 


BE  binding  energy  of  an  adsorbate  relative  to  the  gas  phase 
(eV,  Eq.  (2)) 

AGi  Gibb’s  free  energy  change  of  adsorption  of  BH4-aq  (eV,  Eq. 

(3)) 

K\  equilibrium  constant  for  adsorption  of  the  BH4-aq  (unit¬ 

less,  Eq.  (4)) 

T  temperature  (I<) 

[BH4- ]  dimensionless  concentration  of  BH4-aq  (unitless) 

AG2  Gibb’s  free  energy  change  of  adsorption  of  OH-aq  (eV,  Eq. 

(6)) 

/<i  equilibrium  constant  for  adsorption  of  the  OH-aq  (unit¬ 

less,  Eq.  (5)) 

[OH-  ]  dimensionless  concentration  of  OH-aq  (unitless) 

AG3  Gibb’s  free  energy  change  of  desorption  of  B(OH)3*  (eV, 
Eq.  (8)) 

[B(OH)4- ]  dimensionless  concentration  of  B(OH)4-aq  (unitless) 
I<3  equilibrium  constant  for  desorption  of  the  B(OH)3*  (unit¬ 
less,  Eq.  (9)) 

Ea  potential  dependent  activation  barrier  for  an  elementary 

reaction  (eV,  Eq.  (13)) 

Ea°  activation  barrier  of  an  elementary  step  at  the  equilibrium 

potential  for  H*  desorption  (eV) 

P  symmetry  factor  (unitless) 

/<i  rate  constant  for  dissociation  of  a  B-H  bond  of  BH4*  (s-1 

site-1,  Eq.  (14)) 

Ai  pre-exponential  factor  within  the  rate  constant  for  B-H 
dissociation  of  BH4*  (s-1  site-1 ) 

rn  rate  of  reaction  n(s-1  site-1,  Eq.  (16),  reaction  steps  num¬ 

bered  in  Fig.  1 ) 

0n  fractional  surface  coverage  of  species  n  (unitless,  species 
numbers  in  Fig.  1 ) 

ic  current  in  elementary  units  (electrons  s-1  site-1 ) 

a  area  occupied  by  1  atom  (cm2) 

A  surface  area  of  experimental  electrode  (cm2) 

n  number  of  atoms  per  active  site 

N  total  number  of  atoms  on  experimental  electrode  surface 

S  total  number  of  active  sites  on  experimental  electrode 

surface 

I  current  (As-1) 
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Appendix  B.  Supplementary  data 

Supplementary  data  associated  with  this  article  can  be  found,  in 
the  online  version,  at  doi:  10.1 01 6/j.jpowsour.201 1.07.042. 
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